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f**^ ' In this lecture I discuss jet-shape distributions and describe how from jet evolution one may design Monte Carlo 

' simulations which are used in the analysis of short distance distributions in e+e^-annihilation, lepton-hadron and 

i—i , 1 Introduction 

, A major difRculty is encountered in dealing with perturbative QCD: in computations one uses quarks and 
gluons but only hadrons are observed. Actually, there is no obvious conflict with colour confinement since it is 
impossible to compute distributions with quarks or gluons in initial or final states, they are divergent due to 
ly^ I collinear and infrared singularities. Even, with this general impediment, there are ways to make predictions 
for hadronic distributions. In more than thirty years since the discovery of asymptotic freedom^ an enormous 
amount of information on QCD at short distance has been obtained by perturbative (PT) methods. For 
' instance, collinear (and infrared) divergences are present in DIS structure functions (hadron in initial state) 
lO ' or in e'^e'^ fragmentation functions (hadron in final state). However these divergences factorize so that one 
is able to compute the evolution'^ of these distributions in the hard scale and makes predictions in terms the 
' distribution at a given scale. 

In principle, the only possibility to make absolute predictions (i.e. without involving phenomenological 
Cl^l inputs except for Aqcd) is to compute inclusive short distance distributions such as e+e^ total hadronic cross 
section in which all PT coefficients are finite (however, as I will recall, difhculties come from non-convergence 
of PT expansions). Using inclusive observables one can reach a "complete" description of the emission of 
QCD radiation. To give examples of the amount of information obtained in jet-physics studies I discuss in the 
■ following jet-shape distributions and describe how from jet evolution one may design Monte Carlo simulations. 

2 Jet-shape observables 

In all hard processes (e"''e~-annihilation, lepton-hadron and hadron-hadron collisions) one may introduce a 
large variety of jet-shape observables. Example in e^e~ are thrust T and broadening B 

h h 

with Q the hard scale (total e^e~ center of mass energy). The sum is over all emitted hadrons with p^t 
transverse momentum and "qt rapidity with respect to the thrust axis (which maximizes T). It is clear 
from these examples that different jet-shape observables characterize different aspects of radiation (transverse 
momentum contributes uniformly to B and mostly at large angles to T). In general one finds that jet- 
shape observables {V = r, _B, C, D, ifout, P, Uij • • ■ ) are small for most events. In the hard process under 
consideration, V^— >0 corresponds to the exclusive limit in which the minimum number of hadrons arc emitted 
(two for T or i? in e+e~) so that QCD radiation is characterized by jets around these primary hadrons. 

A quantitative description of jet radiation can be obtained by studying the fully inclusive hadron distri- 
butions for the various jet-shape observables V 

^(>') = E/^e(. -£„,.). (2) 

n \ n / 
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In the PT study of one replaces the sum of hadron momenta with the one of quarks and gluons. Is 

it possible to assume that the PT distributions represent the hadron distributions? A positive answer is 
suggested by the property that in general these inclusive obscrvablcs arc collincar finite (given pi, pj which 
arc coUinear, then V is the same if one replaces them with the single momentum pi + Pi) and infrared finite 
(given Pi, Pj with \pi\ <C \pj\ then V is the same if pi is neglected). Therefore, assuming that hadrons are 
made of collinear or soft partons, V remains the same if one replaces hadrons with partons: hadron flow ~ 
parton flow. In these cases all PT coefficients of T,{V) are finite. One has 

Spt(I^) = a(yX° + 6(FX°+i • • • (3) 

with as = ctsiQ) in MS scheme and uq depending on the number of jets involved (no = 1 for 2-jct observable 
such as T or B). In the exclusive limit y — > all PT coefficient diverge due to collinear and infrared 
singularities. For V ^ these singularities cancel but for small V they leave large logarithmic terms. A 
reliable PT calculation of S(V^) at small V, where it is large, requires the resummation of the most important 
enhanced terms. They are organized as follows 

oo 

lnSpT(F) = 5^{rf„<i"+i + s„<i" + ---}, L = lnV. (4) 

n = l 

In order to control the scale of the logarithms one needs to resum at least dn (double logs, DL) and s„ (single 
logs, SL) terms. To obtain a complete PT prediction for all values of V one has to match in Ept(V) both 
restunmed (4) and the exact (3) expressions. Although resummation procedures are now well established, to 
reach the required SL accuracy one needs very complex tools (Mcllin or/and Fourier transforms, asymptotic es- 
timations of integrals, etc). Recently a numerical program*^ is available that performs automate resummations 
for jet-shape distributions in e+e~, DIS and hadron-hadron collisions. 

Before describing the structure of Feynman diagrams contributing to the enhanced terms in (4) I will 
recall the non-convergence"^ of PT expansion and its relation to the large distance region of confinement. 



2.1 Power corrections 

Various physical facts are at the origin of the non convergence of PT expansions. In general they imply the 
presence of power corrections to PT results. The fact which is phcnomcnologically most important^ is that 
the running coupling is involved at any scale smaller than Q. For example, the average value of V is given by 

{V)= / =^a,{kt)-V{kt/Q), (5) 
Jo '^t 

where the virtual momentum kt in the Feynman diagrams runs into the large distance region (although the 
observable is at short distance). Since the observable is collinear and infrared finite, one has ki^V{kt/Q) ~ Q-^ 
for kt — > 0. Here however the coiipling enters the confinement region. Mathematically this PT difficulty is 
reflected into the fact that, although all PT coefficients in as{Q) are finite, the expansion is non-convergent 
(renormalon singularity). To make a quantitative prediction one has to deal with the large distance region for 
the running coupling. There are various prescriptions for this. 

One prescription^ consists in introducing a non-perturbative (NP) parameter given by the integral of the 
running coupling in the large distance region and expressing (5) as 

(y) = lyfp^ +^^\cv aoit^i) + f2A^aA, a^im) = f" , (6) 

^ [ n=l J -'O 

with /X/ a short distance scale (/U/-independence is ensured by renormalization group) and Cy a known constant 
depending on the observable. Renormalons in the A''-order PT expressions {V)pj, are canceled by the sum 
over the known coefficients Ay. The contribution from the NP parameter Q;o(/i/) is suppressed by inverse 
powers in the hard scale. This power correction is detectable^ even at LEP energies. Similar NP contributions 
are found in the PT study of the distribution T,{V) and also here one needs to introduce the NP parameter 



3 



ao The effect here is in general a power correction "shift" in the argument of the PT result Spt(T^)- A 
consequence of this prescription is that the same NP parameter ao{fii) enters all jet-shape observables and 
then one can study is phenomenological consistence. In general one finds^ that for the various quantities the 
fitted values of ao(/i/) varies within about 20%. 

Higher power corrections are present and maybe important. Here one needs more general prescriptions^ 
with introduction of shape functions to modulate large distance contributions. These prescriptions allows one 
to describe the distributions at low values of V. 



2.2 Structure of PT contributions 

To understand the features of QCD radiation one needs to consider how PT resummation is obtained via 
factorization of (universal) collinear and infrared singularities. As pointed out by Dasgupta and Salam'*, the 
situation is different for global and non-global jet-shape observables. 

Global jet-shape observables. Here the full phase space of emitted hadrons is considered. Examples in 
e+e" are T and B in (1). In these cases the DL and SL contributions in (4) are due to gluon bremsstrahlung 
emission off the primary quark-antiquark. These collinear and/or infrared enhanced contributions factorize 
and are resummed by /mear evolution equations leading to Sudakov form factors. Therefore, after factorization 
of collinear and infrared singularities (including soft gluon coherence) QCD radiations appears as produced 
by "independent" gluon emission. Gluon branching (into two gluons or quark-antiquark pair) enters only in 
reconstructing the running coupling as function of transverse momentum. 

The fact that here the branching component does not contribute (within SL accuracy) can be understood 
as a result of real-virtual cancellations of singularities. Indeed, in the collinear limit, the transverse momentum 
of an emitted gluon is equal to the sum of transverse momenta of its decaying products. Therefore, if one 
measures the total emitted transverse momentum, as in broadening for instance, it is enough to consider the 
contributions of bremsstrahlung gluons (independent emission similar to QED). Further branching does not 
contribute due to unitarity (real-virtual cancellation). 

Non-global jet-shape observables. Here only a part of the phase space of emitted hadrons is considered. 
The best known example is the Sterman- Weinberg distribution^ of energy recorded inside a cone around a jet. 
A simpler example in e^e~ is the distribution in energy recorded outside a cone around the thrust 
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Since the jet region is excluded, there are no collinear singularities to SL accuracy and the resummed PT 
contributions come from large angle soft emission. Here resummation is more complex but informative then 
in the previous cases^'^"'^"'^. Both bremsstrahlung and branching components contribute 






The bremsstrahlung component resums contributions from gluons emitted in the recorded region outside the 
cone. This gives a Sudakov function exponentially decreasing with the SL function 

'■^ dkt N.asjh) 

, kt TT 



(7) 



The branching component resums contributions from gluons emitted inside the jet region. These gluons need 
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to branch in order to generate decaying products entering the recorder region. Here real-virtual cancellation 
is incomplete and virtual enhanced contributions are dominating thus leading to a strong suppression in the 
distribution which asymptotically turns out to be Gaussian in r. 

The described behaviour has been obtained by numerical* and analytical^° studies based on the multi-soft 
gluon emission distributions^^. I describe this calculation which will be used also to introduce Monte Carlo 
simulations for jet physics. 



3 Monte Ccirlo simulation 

I start by describing the derivation of the evolution equation-^° used to compute non-global jet-shape observ- 
ables. One introduces the generating funetional for all multi-soft gluon distributions^^ (known only in the 
planar approximation of large Nc) and shows that it satisfies a branching evolution equation corresponding to 
a Markov process which can be numerically implemented into a Monte Carlo simulation program. To do this 
one needs to include proper cutoff for coUinear and infrared singularities. Prom generated events one computes 
the desired jet-shape distributions (the cutoff contributes only with power corrections). These Monte Carlo 
simulations not only allows one to compute jet-shape distributions, but provides a "complete" description of 
the hard process with soft gluons emitted in any angular region (including large angles). The next crucial step 
(apart for including quarks and non soft contributions) consists in how to go from parton to hadron emission. 



3.1 Generating functional (soft and planar limit) 

The starting point is the amplitude for the emission of n-soft gluon ki, - ■ ■ ,k„ off a primary qq pair of 
momentum p, p. It is represented as a sum of Chan-Paton factors and the colour ordered amplitudes coefficients 

Mnippqi ■■■qn)= i^""' ■ • • }/3^ 

perm. 

Prom the factorization of the softest emitted gluon q-^ 

Mr,{-..£ni'...)=g,M^_,{-..£i'-..).(j^-j^), (9) 

\[qiqn) {qi'qnjj 

one deduces a recurrence relation and computes all colour amplitudes in strong energy ordering (leading order 
in soft limit). Summing over colour and polarization indices, the distribution is given, to leading TVc-order, by 



perm. 



where {ij) = l — cos6ij and the emission is off a general o6-dipole. Similar approximations give the multi-soft 

gluon distributions in pure Yang-Mills theory. It is interesting that this multi-soft gluon distribution coincides 
with the square of the exact MHV colour amplitude discussed by Parke- Taylor^^. 

To study arbitrary jet-shape distributions one introduces a source function u{q) for each soft gluon 

^,,{E, u)=Y,-,l— j{<<ii) = E /n(— 4^ 4 • •••'")' (11) 

with = Ncas/n and E = Q. This functional summarizes the full information for the soft gluon emission. It 
involves only real emission distribution so that it is infrared and coUinear divergent. Virtual corrections will 

be included later at the same accuracy in the soft limit. 

The evolution equation is obtained by using the factorization structure of multi-soft gluon distribution 

Wab{l,--- ,n) = j^^^Wai{l,--- ,^-l)-W^^6(^+l,--- ,n), (12) 
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with one of the soft gluons. Taking qi as the hardest (soft) gluon one obtains 

EdE^ab = J ^^asWabiq) [u{q) Y^aq ■ Yqb - '^ab] , Wabiq) = (-^^"^^^^ (1^) 

The negative terms in the integrand corresponds to the virtual corrections (via Cauchy integration) . They are 
included in the scheme in which, for the fully inclusive case of u{q) = 1, they completely cancel against the 
real contribution {T,{E, 1) = 1). Both the real emission branching (first term in the integrand) and the virtual 
correction are coUinear and infrared singular. For inclusive observables, (i.e. for suitable sources u{q)) these 
singularities cancel. 



3.2 Monte Carlo simulation, soft gluons at large angles 

The evolution equation (13) can be formulated as a Markov process and then numerically solved by Monte 
Carlo simulations®. Here the basic system is the a6-dipole which, emitting a soft gluon q, branches into the 
two dipole aq and qb. To construct a Mante Carlo process" one rewrites (13) by splitting real and virtual 
corrections. To do so it is necessary to introduce a cutoff Qo in transverse momentum (the argument of as). 
The Sudakov form factor 

rip r HQ 2E'^ 
\nS^,{E) = - / ^ / -^asW^M ■ e{q,^,-Q,) , ql, = , (14) 

JQo J Wab[q) 

is the solution of (13) in which the real emission piece in the integrand is neglected. Here qtab is the transverse 
momentum of q with respect to the a6-dipole. Then one introduces the probability for dipole branching: 
{ab)^{aq){qb), E^E^ 

dV{E„ O,) = { ^ -g^ } { ^aswM] ■ 0{qtab-Qo) • (15) 

Here the first factor selects Eg and the second the direction of the emitted soft gluon. The branching 
does not take place if the momentum qtab is below the cutoff. Each emitted dipole then undergoes successive 
branchings with decreasing energy until no further branching is permitted by the cutoff. Each Monte Carlo 
run generates an "event" with soft gluon emitted above the cutoff. Unitarity of probabilities ensures, in fully 
inclusive distribution (M(g) = l), the complete cancellation of real-virtual contributions. Computing jet-shape 
observables one has cancellations of collinear and soft singularities with residual In ^^-contributions and power 
corrections in Q()/{VE). The accuracy reached in the calculation of jet-shape observable is based on the fact 
that here one uses (planar) multi-gluon emission distributions. This implies that for jet-shape distributions 
one has included both DL and SL of soft origin. SL terms of collinear origin are missing so that here the 
g^gg splitting function 

Pg^M=Nc{^^ + ^^+z{l-z)-2^, (16) 

includes only the infrared singular pieces for 2;— s-0 or 2; ^ 1. Together with the non-soft pieces of the splitting 
function also quark branching channels are missed here. 

An additional crucial missing clement for a realistic jet-emission simulation is the fact that no dipole 
momentum recoil is here taken into account. However this branching formulation correctly accounts for soft 
emission also at large angle which are contributions missed in the present Monte Carlo simulation. It would 
be then important to include recoil in (15). 



"A Monte Ccirlo program based on dipole branching similar to what is described here has been constructed by L. Lonnblad 
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3.3 Im,proved Monte Carlo simulation 

The most accurate available Monte Carlo simulations^^' are based on a branching algorithms'^ '^^ which 
resum (for instance for jet-shape distributions) DL and some of the most important SL terms. In particular 
soft gluon coherence is included (angular ordering). Continuous upgradings of Monte Carlo codes are underway 
which account for new theoretical, phcnomcnological and experimental results. I list here** some of the major 
features and recent (or future) developments. 

Full collinear singularity structures are included. Branching is in general formulated as successive parton 
emission'"''''^ so that recoil and splitting functions for all parton branching (g gg, g qq. q qg) arc 
naturally included. A major development in the structure of branching is expected including resummation of 
soft radiation at large angles. Partial results are however available^^. 

Parton branching discussed in this lecture involves only final state emission. Similar branching holds for 
initial state radiation needed in hadron-lepton or hadron-hadron collisions. However in the small- .t region (i.e. 

<C s) there are peculiar differences. Here soft gluons are both emitted and exchanged partons. For soft 
exchanged gluons angular ordering requires^^ additional virtual corrections of non-Sudakov type (connected to 
gluon "Reggeization" ) . This brings one into the domain of "small- a; physics" in which an accurate branching 
algorithms has still to be found. Partial results are however available^°. 

Using factorization structure of collinear and infrared singularities the same branching structure holds for 
all hard processes. It is then possible to construct a single Monte Carlo code for e+e~ annihilation, lepton- 
hadron DIS and hadron-hadron (large Ex)- Moreover, such a factorization structure allows one to include in 
the QCD Monte Carlo code also non-QCD processes (Electro-weak, beyond the standard model, gravity,...). 
Such Monte Carlo codes are then very useful instruments for quantitative study of "new-physics" scenarios. 

Major recent developments in the Monte Carlo codes are systematic attempts^' to account for NLO and 
NNLO exact results including heavy flavour processes. 

3.4 From partons to hadrons 

The above description of the Monte Carlo code refers to the generation of events with emission of partons 
(possibly together with non-QCD particles). As noted before, due to the presence of collinear and infrared 
singularities these emission processes require a cutoff Qo- The question is then how to go from partons to 
hadrons and how this "affects and distorts" the QCD radiation. Hadronization in the Monte Carlo code is 
based on the property of preconfinement^^: after successive branchings, partons are emitted in clusters of 
colour singlets of mass of order Qq. It is then natural to convert these colour singlet clusters into hadrons 
without "affecting and distorting" the QCD radiation within a scale of order Qo- Preconfinement may be 
basis for the property that parton flow ~ hadron flow 

The basis for preconfinement is again the collinear and infrared structure of QCD. It can be explained as 
follows. In the large Nc limit, one may follow the colour line of partons in the successive branchings. The 
colour line of an emitted quark (or quark part of a gluon) ends into the colour line of an emitted antiqiiark 
(or antiquark part of a gluon). Due to this colour connection, they form a colour singlet. If two partons are 
colour connected, no emission takes place along the colour line which connects them. So virtual corrections are 
dominating and the distribution in the mass of the two colour connected partons is suppressed by a Sudakov 
form factor. As a result the mass of the two parton system is of order Qo- 

4 Final considerations 

Since 1973 exact and resummed PT calculations produced enormous inside on QCD radiation and its phc- 
nomcnological "evidence" . A clear sign of this success is given by the Monte Carlo codes which summarize 



''Of course the following list does not really account for all the work done in the field. It reflects my personal view of the important 
points. 
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many QCD results and are used for the analysis of data in all hard processes (within known theoretical accu- 
racy). Their phenomenological success is also an indication that it is "reasonable" to assume that parton flow 
~ hadron flow. 

In PT studies it is possible to circumvent the difficulties of colour confinement and hadronization via 
factorization properties or "reasonable" phenomenological assumption. However this is not satisfactory on the 
theoretical point of view. There are indications that the problem of colour confinement could be approached 

by a dual formulation of QCD in terms of extended objects. Lattice QCD calculations indicatc^^ that QCD 
vacuum is populated by extended low dimensionality objects responsible for confinements. String theory 
provides in principle a basis for this study. Here there are attempts^^ to develop new ideas for PT calculations. 
Also Regge behaviour of scattering amplitudes is studied^^ and the language has some similarity with the 
"Poneron" prc-QCD topological cxpansion^^. However the natural questions arc how NP formulations could 
account for the enormous PT "evidence" and how they could inspire modeling hadronization. 
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